How long-range interactions tune the damping in compact stars 
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Long-range interactions lead to non- Fermi liquid effects in dense matter. We show that, in contrast 
to other material properties, their effect on the bulk viscosity of quark matter is significant since 
they shift its resonant maximum and can thereby change the viscosity by many orders of magnitude. 
This is of importance for the damping of oscillations of compact stars, like in particular unstable 
r-modes, and the quest to detect signatures of deconfined matter in astrophysical observations. We 
find that, in contrast to neutron stars with standard damping mechanisms, compact stars that 
contain ungapped quark matter are consistent with the observed data on low mass x-ray binaries. 
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Quantum chromodynamics (QCD) features strong, 
long-range gauge interactions that lead to color con- 
finement. At high density these interactions are partly 
screened by the medium and can lead to deconfine- 
ment and various hypothetical forms of quark matter 
[T]. In the color-flavor locked (CFL) phase at asymp- 
totically large density all quark modes are gapped, but 
many of the possible phases at moderate densities that 
might be reached in compact stars feature ungapped 
quark excitations. Yet, in-medium effects cannot entirely 
screen the gluonic interaction, but their magnetic part is 
only Landau-damped These long-range correlations 
strongly modify the low-energy fermionic excitations of 
deconfined dense matter and lead to non-Fermi liquid 
(NFL) behavior [31 S] of various material properties |51- 
[S] . Analogous NFL effects due to long-range interactions 
are studied extensively in condensed matter systems [5] 
and with dual holographic descriptions [lOj. Whereas 
in such systems NFL effects are generally only found in 
controlled experiments at low temperatures, we discuss 
here a case where they have important observable conse- 
quences outside of a physical laboratory and in an envi- 
ronment hotter than the surface of the sun. 

Compact stars are the only known objects that are 
dense enough that they could contain deconfined quark 
matter, and their mechanical oscillation modes are the 
only known way to directly probe the dense matter in 
their interior. The connection is established via the 
damping of these modes which is locally described by 
the viscosities of the particular microscopic form of mat- 
ter. Whereas the shear viscosity is generally very similar 
in hadronic and quark matter, since both arise from un- 
screened plasmon-enhanced scattering of relativistic par- 
ticles, namely leptons jTT] and ungapped quarks [S], re- 
spectively, the bulk viscosity can be decisively different 
for various forms of matter |12ffT7] . In dense matter bulk 
viscosity arises generally from slow flavor-changing weak 
interactions. The bulk viscosity is a resonant mechanism 
that features a sharp maximum when the time scale of 
the relevant weak process matches the one of a driving 
density oscillation. Although the bulk viscosity is in- 
duced by weak processes, their rate depends on strong 
interactions since they can open phase space for weak 



transitions. Correspondingly low energy strong interac- 
tions control the resonance, analogous to the dial of a 
radio tuner, and even moderate changes can alter the 
viscosity at a given temperature and frequency by many 
orders of magnitude. 

In this work we study the influence of NFL interac- 
tions on the bulk viscosity of dense quark matter. The 
bulk viscosity stems from a lag between a driving density 
oscillation and the chemical equilibration of the system 
mediated by slow weak interactions. Here we study the 
damping of small amplitude mechanical oscillations in 
the subthermal case where the displacement of the chem- 
ical potentials /^a from their equilibrium values satisfles 
/iA ^ T. The suprathermal regime [131 [TB] and the in- 
fluence of strong interactions on the saturation of unsta- 
ble modes jTH] will be discussed elsewhere. Whereas in 
dense hadronic matter chemical equilibration is caused 
by semi-leptonic Urea processes, in dense quark matter 
strangeness-changing non-leptonic processes d+u ^ u+s 
dominate [13 . In this case the displacement fi/^ — fid—fJ-s 
is finite when density oscillations drive the system out of 
chemical equilibrium and causes a net non-leptonic rate 
F^^-* = Tii^g — Tg^^ that re-equilbrates the system. To 
compute this rate we follow the analysis performed in 
and extend it to the interacting case. The computa- 
tion requires the dispersion relations of low energy quarks 
which is modified by self-energy corrections from quan- 
tum fluctuations involving unscreened magnetic gluons. 
The form of the dispersion relation including such NFL 
effects is close to the Fermi surface given by [3J [TJ HU] 
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where in the weak coupling, hard dense loop |21| approx- 
imation the parameters are given by [6l [7] 
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The strong correlations signiflcantly increase the density 
of states near the Fermi surface where the group veloc- 
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ity of the corresponding excitations gradually deviates 
from its Fermi liquid value vp and eventually vanishes in 
the low energy limit. This effect arises from a kinematic 
low energy enhancement and is even present in the weak 
coupling regime. Yet, the parametric form of the NFL 
dispersion relation has been shown to hold beyond per- 
turbation theory [2] and via a manifest power counting 
analysis it has been pointed out that this result holds 
within an effective theory in the low temperature limit 
even at strong coupling fl| . Aside from the strong generic 
increase of the density of states close to the Fermi sea 
due long-range interactions, the other strong interaction 
corrections to the Fermi liquid parameters are not log- 
arithmically enhanced and due to our ignorance on the 
values of these parameters in the strong coupling regime 
relevant for compact stars we neglect these moderate cor- 
rections here. As discussed in ^ there are likewise no 
logarithmic corrections to the weak interaction vertices. 

In the subthermal limit /iA ^ T <C /i and neglecting 
the subleading quark masses mj ^ /i we find for the 
non-lcptonic rate 
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This result differs from the non-interacting expression 
given in |13l [TO] by the strong interaction corrections in 
terms of NFL factors A = l-|-tTlog(A/T) for each of the 
involved four quarks. These factors become large in the 
ultradegenerate regime T A = 0(/x) realized in com- 
pact stars. In eq. (|3| additional terms suppressed in 
powers of cr/A which could not be evaluated in closed 
form have been neglected, but a numerical integration 
shows that their combined contribution is at the relevant 
temperatures less than 5%. 

For comparison we will also study the viscosity due to 
Urea processes including NFL interactions. The direct 
Urea rate is obtained as a generalization of the Fermi 
liquid expression \ri\ [^ . In case of the dominant d- 
quark Urea processes d ^ u + + Dg, u + d + Ve 

where /iA — fJ-d — fJ-u — fJ-e we find for the NFL dispersion 
relation eq. ([ij 
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As observed previously in the study of the neutrino emis- 
sivity [7 , this expression exhibits a similar but weaker 
NFL enhancement since only two quarks are present in 
the weak process. The same holds for the Cabbibo sup- 
pressed s-Urca channel. However, these Urea rates are 
strongly suppressed in r//i compared to the non-leptonic 
rate eq. (|3]). 

The bulk viscosity is a measure for the local dissipa- 
tion caused by a (harmonic) oscillation of the conserved 
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Table I: Strong interaetion parameters, defined in eq. 
describing the response of the particular form of dense mat- 
ter. For comparison we also show the values for interacting 
hadronic matter in terms of the proton fraction x and the sym- 
metry energy 5*. The corresponding weak rate for hadronic 
Urea processes is given in |12) . 



baryon density n(f, t) = n{r) + An,(f^ sin(a;t) and can be 
defined in terms of the energy dissipation by 
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Its detailed derivation is for instance discussed in [16J. 
Generally the non-leptonic and the Urea channels mix 
and require a coupled-channel analysis |22l [25] , but due 
to the smallness of the Urea rate this effect is negligi- 
ble and we evaluate them here individually. In the sub- 
thermal limit, denoted by the superscript "<", the bulk 
viscosity has the general resonant form [16] 
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where the resonant maximum (^^ 

pendent of the weak rate 7 = F^^^ /ma, but the latter de- 
termines at which temperature this maximum is reached. 
Here the strong susceptibilities B and C determine the 
response of the medium and are given by 



C = 



dn 



B = 



h dx 



(5) 



where x = n/^/n is the fraction of the particle density 
corresponding to /iA that is driven out of chemical equi- 
librium. The susceptibilities are given for the different 
processes in table [Tj and in particular they do not involve 
logarithmic corrections. This can be seen from the pres- 
sure at finite temperature and density which has been 
computed in weak coupling beyond leading order in f24]. 
This result shows that all logarithmic temperature cor- 
rections are further suppressed in T/ fi. This is generally 
expected since in the pressure and the susceptibilities the 
entire Fermi sea contributes and the contribution from 
the thin shell of width 0(T) around the Fermi surface is 
subleading. 

Let us first consider the expected size of the NFL 
corrections. From our knowledge of the running of the 
strong coupling |25| at scales 0(/i), quark matter at neu- 
tron star densities is, analogous to the strongly correlated 
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Figure 1: The influence of NFL efi'ects due to long-ranged 
magnetic gauge interactions on the bulk viscosity of dense 
quark matter. All curves are given for a baryon density 
n = 3no and a frequency u = 1 kHz. Shown is the dominant 
contribution due to non-leptonic processes for three values of 
the strong coupling constant {as=0 dotted, as = 1 solid and 
Os — 3 dashed), as well as the subleading contribution from 
d-quark Urea reactions. For comparison the bulk viscosity 
due to modified Urea interactions in dense APR neutron star 
matter, studied in |16| . is shown as well. 



Figure 2: Impact of the NFL corrections to the bulk viscosity 
on the r-mode instability region for a 1.4 A/q strange quark 
star. The latter is shown for two values of the strong cou- 
pling constant {as —0 dotted and as = l solid). For compar- 
ison the instability region of a 1.4 M© neutron star is shown 
by the dot-dashed curve. For details on the star models see 
[26| . The theoretical curves are compared to LMXB pulsar 
data compiled in [27|, where either temperature estimates are 
known (points) or only upper bounds were possible (points 
with arrows). 



fluid observed at RHIC, expected to be strongly coupled 
as = 0(1). Using the perturbative expressions eq. (2) 
as a guide, the relevant range of the parameters eq. (2) 
in compact stars should be roughly 0.05 ^ cr < 0.3 and 
100 < A/T < 10*^ so that the enhancement factor A"* 
can take values in the range 2 < X'^ < 1000. Whereas 
the logarithmic form eq. ([T|| is generic, the perturbative 
relations between the parameters a and A and the cou- 
pling do not need to hold, but due to our ignorance of 
the actual dependence in the strong-coupling regime we 
extrapolate them below. 

Our results for the bulk viscosity are shown in fig. [T] 
for different values of the strong coupling, using an inter- 
mediate value for the strange quark mass m^. = 150 MeV. 
As can be seen the increase of the non-leptonic rate due to 
long-range gluonic interactions shifts the resonant maxi- 
mum of the bulk viscosity to significantly lower tempera- 
tures and thereby strongly increases the damping of cold 
quark matter. As observed before in |22l [53] the Urea 
contribution is only relevant at small frequencies and 
temperatures T > 10^ K. Due to the moderate change 
of the Urea rate this conclusion is unaltered by NFL ef- 
fects. Compared to the corresponding bulk viscosity in 
hadronic matter the damping of interacting quark matter 



is more than ten orders of magnitude larger at tempera- 
tures T < 10^ K present in cores of old compact stars. 

As an example for the importance of the considered 
strong interaction corrections to the bulk viscosity of 
quark matter we study the instability of r- modes p8l [29] 
of strange quark stars |3U]. R-modes are global oscil- 
lation modes of stars rotating at a finite angular veloc- 
ity n, that are unstable to the emission of gravitational 
waves within a characteristic instability region in the T- 
il-plane. At the boundary of the instability region the 
gravitational instability is precisely balanced by viscous 
damping whereas for frequencies above the r-mode grows 
exponentially, for more details see |26j . If r-modes are 
unstable they strongly radiate gravitational waves and 
would quickly spin down a star [31]. Therefore, if the 
r-mode growth is stopped at values predicted by known 
saturation mechanisms |18l [32] , observed low mass x-ray 
binaries (LMXBs) , that are slowly spun up over long time 
periods by mass transfer from a companion, should not 
lie in the r-mode instability region. Yet, as can be seen 
in fig. [2] most observed LMXBs lie clearly in the in- 
stability region of a 1.4 Mq neutron star with standard 
damping mechanisms, i.e. shear viscosity due to leptons 
[TT] and bulk viscosity due to modified Urea reactions 
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|12| . The situation is even worse for heavier stars and 
this conclusion cannot be changed by the involved uncer- 
tainties due to a remarkable insensitivity to the underly- 
ing parameters [25]. The enhanced dissipation of phases 
with ungapped quarks strongly increases the damping at 
low temperatures and leads to a stability window where 
r- modes are stable up to large frequencies. Neverthe- 
less, neglecting interactions in the bulk viscosity many 
observed LMXBs are still within the instability region 
of a 1.4 Mq strange star, as shown by the dotted curve. 
However, taking into account long-range gauge interac- 
tions in the strong coupling case as — I all observed pul- 
sars are consistent with lying outside of the instability 
region. Therefore, the enhanced damping of ungapped 
quark matter can explain the large spin rates of LMXBs. 

The results found here for strange stars immediately 
generalize to hybrid stars, where a quark matter core 
is enclosed in a hadronic outer part. Since the compo- 
nent with the stronger viscosity dominates the damping 
the instability region of a hybrid star is roughly given in 
each case by the upper segments of the neutron star and 
strange star curves • Further our results apply to cer- 
tain color superconducting phases that feature ungapped 
quark species and unscreened gluonic modes that can in- 
duce a NFL enhancement, see [7]. However, in hyper- 
onic stars [33 where similar non-leptonic processes exist, 
there is no NFL enhancement since long range color in- 
teractions are absent. 

It is striking that there are no data points for pulsars 
spinning with frequencies of several hundred Hertz that 
at the same time have core temperatures < 10^ K. A sta- 
bility window, as found here for ungapped quark matter, 
provides a robust explanation for this observation since 
the strong dissipative heating [34' due to r-modes in the 
instability region prevents stars from entering it. Clearly, 
our results involve considerable uncertainties on the de- 
tailed strong coupling values of the various low energy pa- 
rameters entering the bulk viscosity. Nevertheless, within 
the uncertainties interacting quark matter can provide 
a consistent explanation for the the astrophysical data, 
whereas hadronic matter with standard damping mech- 
anisms can not [26J. Potential enhanced damping mech- 
anisms in hadronic stars related to the complicated neu- 
tron star crust have been proposed, like boundary layer 
rubbing or mutual friction in hadron superfluids, see e.g. 
[27]. It requires detailed further study to decide if these 
mechanisms can likewise explain the astrophysical data, 
or if the data might eventually provide a robust signa- 
ture for deconfined quark matter. In this endeavor the 
connection of the dynamic star evolution to pulsar tim- 
ing data [34_ should provide important information to 
distinguish between these different scenarios. 
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